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We discuss the treatment of intrinsic transverse momenta in high energy scattering processes. 
Within the field theoretical framework of QCD the process is described in terms of correlators 
containing quark and gluon fields. The correlators, parameterized in terms of distribution and 
fragmentation functions, contain matrix elements of nonlocal field configurations requiring a careful 
treatment to assure color gauge invariance. It leads to nontrivial gauge links connecting the parton 
fields. For the transverse momentum dependent correlators the gauge links give rise to time reversal 
odd phenomena, showing up as single spin and azimuthal asymmetries. The gauge links, arising from 
multi-gluon initial and final state interactions, depend on the color flow in the process, challenging 
universality. 
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I. INTRODUCTION 

The basic degrees of freedom that feel the strong interactions, quarks and gluons, are confined into hadrons, strongly 
interacting particles. Considering the nucleons (light hadrons), the characteristic energy and distance scales are given 
by the nucleon mass Mn, or taking into account the color degrees of freedom one may prefer a scale Mn/N c ~ 300 
MeV. We refer to this as 6(M) or 0{Q Q ) if we consider high-energy processes. Such processes are characterised by 
hard kinematical variables that are of order Q with Q 2 ^> Afjy. Depending on details, the high-energy scale Q can be 
the CM energy, Q ~ y/s or it can be a measure of the exchanged momentum. 




FIG. 1: Schematic illustration of the contribution of a hard subprocess, parton (pi) + parton (P2) — * parton (fci) + parton 
(fe), to the (2-particle inclusive) scattering process hadron (Pi) + hadron (P2) — > hadron (Ki) + hadron (K2) + X, at the 
level of the amplitude. The process being hard implies for the hadronic momenta Pi ■ P2 ~ Pi ■ K\ ~ Q 2 , etc. 



The basic framework for the strong interactions is Quantum Chromodynamics (QCD). Hadrons, however, do not 
correspond to free particle states created via the quark and gluon operators in QCD. The situation thus differs from 
that of QED with physical electrons and photons. In the latter case one knows how in the calculation of an S-matrix 
element contraction of annihilation and creation operator in the field and particle state lead to the spinor wave 
function. For positive times £° = t, one has 

<0|Vi(Ob> = b\p)\0) = (0\M0)\P) e~' ip -« - «i(p) e" 4 ^, (1) 
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with p° = E p = ^Jp 2 



Such a matrix element is 'untruncated' as seen e.g. from 



(o\Mo\p)m=m 



d 4 k 



,-tfe-e 



fc 2 - 



ze 2to 



(27r) 3 2£;p(5 3 (fc-p). 



(2) 



In a process involving a composite hadronic state \P), contractions with one or several of the quark and gluon operators 
may be involved, leading to nonzero matrix elements for a quark between the hadron state and a remainder, but also 
to nonzero matrix elements involving multi-parton field combinations, 

(x\MO\P), WW(£)|p>,... . 

For a particular hadron and a parton field combination, one may collect those operators that involve hadron \P) into 
(distribution) correlators 



®ij(p;P) = £ / 

X J 



d*P x 



PWA0)\X) Wi(0)|P) S\p + Px-P) 



1 



(2tt)< 



(2tt) 3 2E x 

d 4 C (P|^.(0)^(£)|P), 



or correlators involving matrix elements of the form 
1 



*ft(p,Pi;P) = 



(2tt) 8 



(3) 



(4) 



pictorially, 




*(p;P,S) 




P-P, 



^(p,p- Pl ;P,S) 



We will not attempt to calculate these, but leave them as the soft parts, requiring nonperturbative QCD methods 
to calculate them. In particular, although being 'untruncated' in the quark legs, they will no longer exhibit poles 
corresponding to free quarks. These are fully unintegrated parton correlators for initial state hadrons, in general quite 
problematic quantities. For example, they are by themselves not even color gauge-invariant, an issue to be discussed 
below. When more hadrons are involved, one could consider two-hadron correlators, involving two-hadron states (or 
correlators involving hadronic states in initial and final state), etc. If the hadrons are well-separated in momentum 
phase-space with Pi ■ Pj ~ Q 2 , one expects on dimensional grounds that incoherent contributions are suppressed by 
l/(Pj — Pj) 2 ~ 1/Q 2 and one can (at least naively) factorize using forward correlators for single hadrons, connected 
by a hard partonic subprocess. Such a separation in momentum space requires a hard inclusive scattering process 
(Q 2 ~ s). The inclusive character is needed to assure that partons originate from one hadron, leaving a (target) 
jet. In turn, final state partons decay into a jet, in which we also consider a single identified hadron, which can 
straightforwardly be extended to a multi-particle, e.g. two-pion, state. For the fragmentation process of a parton 
(with momentum fc) into hadrons (with momentum P^) we combine the decay matrix elements in the (fragmentation) 
correlator, for quarks 



A„(fc,P ft ) 



£ 

X 



1 



d^ e lk < (Q\^ i (O\P h ,X)(P h ,X\^ j (0)\0) 



(27T) 



d 4 £ e lk < (0|V^Ka^,(0)|0>, 



t„ rr. 



(5) 



where an averaging over color indices is implicit. Pictorially we have 



A(k;P h ,S h ) 



- k 
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In particular, we note that in fragmentation correlators one no longer deals with plane wave hadronic states, but with 
out-states \Ph, X). In all of the hadronic states mentioned before one can also consider polarized hadronic states. The 
spin of quarks is contained in Dirac structure and that of gluons in the Lorentz structure of correlators. 

The basic idea in the diagrammatic approach is to realize that the correlator involves both hadronic states and 
quark and gluon operators. The correlators can be studied independent from the hard process, provided we have 
dealt with the issue of color gauge invariancc. The correlator is the Fourier transform in the space-time arguments 
of the quark and gluon fields. In the correlators, all momenta of hadrons and quarks and gluons (partons) inside the 
hadrons are soft which means that p 2 ~ p ■ P <~ P 2 = <C Q 2 ~ s. The off-shellness being of hadronic order 
implies that in the hard process partons are in essence on-shell. Consistency of this may be checked by using QCD 
interactions to give partons a large off-shellness of &{Q) and check the behavior as a function of the momenta. In 
these considerations one must also realize that beyond tree- level one has to distinguish bare and renormalized fields. 



II. 



COLLINEAR AND TRANSVERSE MOMENTUM DEPENDENT CORRELATORS 



In a hard process, the parton fields that appear in the different correlators correspond to partons in the subprocess 
for which the momenta satisfy pi ■ pj <~ Q 2 . In the study of a particular correlator it implies the presence of a 'hard' 
environment. To connect the correlator to the hard part of the process, it is useful to introduce for each correlator 
with hadron momentum P, a null-vector n, such that P ■ n ~ Q. Using this relation, n would be dimensionless. It is 
actually more convenient to replace n/{P-n) by a dimensionful null-vector n ~ such that P-n = 1. The vectors 
P and n can be used to keep track of the importance of various terms in the correlators and in the components of 
momentum and spin vectors. If one prefers a dimensionless vector, one must make a choice P-n ~ Q. In that case all 
further appearances of n in this section should simply be replaced by n/(P ■ n). The n- vector will acquire a meaning 
in the explicit applications or play an intermediary role. At leading order, it will turn out that the precise form of n 
doesn't matter, but at sublcading (1/Q) order one needs to be careful. 

For parton momenta we write the Sudakov decomposition 



p = x P + p T + (p ■ P - x M 2 ) n, 



(6) 



M 2 /Q. We have the exact relations p ■ p T 



p 2 T 



{p-xPf. The 



where the term x P ~ Q, p T ~ M and an 
momentum fraction x — p ■ n is ^(1). 

In a hard process, the different importance of the various components allows up to specific orders in 1/Q, an 
integrations over some components of the parton momenta. The fact that the main contribution in <I>(p; P) is assumed 
to come from regions where p ■ P < M 2 , whereas the momenta have characteristic scale Q, allows performing the o- 
intcgration up to M 2 /Q 2 contributions (and possible contributions from non-intcgrablc tails). The resulting transverse 
momentum dependent (TMD) correlators are light-front correlators, 



®ij(x,p T ;ri) 



J d(p ■ P) 



d{i-P)d 2 Z T lp< 
(2tt)3 



<P|^.(0)^(OI^) 



(7) 



LF 



where we have suppressed the dependence on hadron momentum P. The subscript LF refers to light-front, implying 
£ • n = 0. The light-cone correlators are the correlators containing the parton distribution functions depending only 
on the light-cone momentum fraction x, 



J d{p ■ P) 



d 2 p T $ij(p;P) 



(27T) 



(8) 



LC 



where the subscript LC refers to light-cone, implying £ • n = (, T = 0. This integration is generally allowed (again up 
to M 2 /Q 2 contributions and contributions coming from tails, e.g. logarithmic corrections from 1/p 2 tails) if we are 
interested in hard processes, in which only hard scales (large invariants ~ Q 2 or ratios thereof, angles, rapidities) are 
measured. If one considers hadronic scale observables (correlations or transverse momenta in jets, slightly off-collincar 
configurations) one will need the TMD correlators. 

The correlators encompass the information on the soft parts. They depend on the hadron and (contained) quark 
momenta P and p (and spin vectors). The structure of the correlator is reproduced from these momenta incorporating 
the required Dirac and Lorentz structure. Clearly, it is advantageous to maximize the number of components along 
the momentum (collinear). For the soft scalar objects this means maximizing the number of contractions with n. This 
leads for nonlocal operators to the dominance of the twist- 2 operators 



and 



pna pn/3 ^ 



(9) 
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(the latter with transverse indices a and 0). Twist in this case is just equal to the canonical dimension of the operator 
combination (remember that dim(n) = -1). 

Of course the appearance of the field strength tensor rather than the gauge field is a requirement of gauge invariance. 
Besides the field tensor, we need the inclusion of gauge links 



U§\ j = 9 exp {-% jf* d( V ■ P) n ■ A(r?)J , 



(10) 



connecting colored parton fields. In the case of the collinear correlators, the gauge links can be built from the 
gauge fields A + — A n — n ■ A, giving a link along the light-cone (£ + = n ■ £ — £ T — 0). The color gauge- invariant 
light-cone correlators for quarks and gluons are 



LC 



T^(x;n) = J ^^—L e*** (P|Tr (F^(0) F™(£) t/gj \P) 
Using the Taylor expansion of the color gauge-invariant nonlocal operators, 



(11) 
(12) 



LC 



,/>t 



^(o)^ 1 ...iD^ V(0) 



one recovers the irreducible set of symmetric traceless local operators relevant in the Operator Product Expansion 
(OPE) approach to describe §%j{x) and T a p(x), namely 



quarks ij 



V' 7 (0)7 {M1 iD» 2 . ..iD^ } ipi(0) - traces, 



0&£ aP = -F^(0)iD^...iD^F^l(0)- traces, 



in which the spin n represents the number of symmetrized indices. Subtracting traces is needed to have an irreducible 
set. The TMD light-front correlators 



^ ] {x,p T ;n) = 
r a/3 \x,p T ;n) = 



d(£ ■ P) d 2 ^ 
(2tt) 3 
d(^P)cf^ 
(2tt) 3 



[P\ 



\^) u m^)\P) 



LF 



(P|Tr (F"^(0) U^f F\(0 Ulffi) \P) 



(13) 
(14) 



LF 



involve a more complex link structure, leading to a path dependence in the definitions (indicated by the arguments 
C and C). This arises because of the (necessary) transverse piece(s) in the gauge link. The simplest possibilities for 
the links in the case of quark and gluon correlators are shown in Fig. [5] [l|. 



III. THE OBSERVABLES 



An important aspect of incorporating intrinsic transverse momenta is the possibility to access them in experiments. 
Consider the subprocess J*(q) + q(j>) — > q(k) which dominates at leading order the inclusive deep inelastic scattering 
(DIS) process 7* (q) + N(P) — > X. In collinear approximation (p w xP) one finds for the momentum fraction x the 
well-known relation x = xb = Q 2 /2P ■ q, i.e. the fraction is identified with the Bjorken scaling variable. For the 
semi-inclusive deep inelastic scattering (SIDIS) process j*(q) + N(P) — > h(Ph) + X with the same subprocess one 
has in collinear approximation (k ss Ph/z) the relation z = Zh — P ■ Phi P ■ q- The non-collinearity in this process is 
determined by 

q T = q + x B P-—P h , (15) 

Zh 

which is taken zero in collinear approximation. First of all, we note that q T is experimentally measurable being a 
difference of vectors of &{Q). This difference is meaningful at ff(M) or higher because mass corrections, coming from 
the identification of x with xb and z with Zh, appear only at ff{l/Q 2 ). The vector q T is the transverse momentum 
of q in a frame in which P and Ph are chosen parallel or (experimentally more useful) it is the transverse momentum 
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FIG. 2: Simplest structures (without loops) for gauge links and operators in quark correlators (a)-(b) and gluon 
correlators (c)-(f). 



—Phx/zh in a frame in which q and P are chosen parallel. We write Q\ — —q T - When Q T ~ ^(M) one easily sees 
that the Sudakov expansions for the quark momenta, p s» + p T and few ^ Ph + k T , imply that q T — k T ~ p T and 
the ^-dependence is attributed to the (convoluted effect of the) intrinsic transverse momenta in the fragmentation 
and distribution correlators. When Q T ~ @{Q) a collinear description involving a subprocess with one additional 
parton radiated off is needed, but for consistency one also wants a match with the TMD description Q ■ 

Not only in electroweak processes like SIDIS or the Drell-Yan process transverse momenta can be accessed. This 
can also be achieved for hadron-hadron scattering. Also here the identification of the transverse momentum is only 
possible together with the identification of the hard subprocess (like pictured in Fig. |T|) . We define 

q T = ^ Ki + ± K 2 - Xi Pi - x 2 P 2 = Pit + P2t - k\ T - k 2T , (16) 

a relation valid up to 0{M). The momenta involved to find q T are in principle all 6 iff) and using at leading order 
9 T «0 yields relations for the momentum fractions in terms of the external hadron momenta (up to 1/Q 2 corrections). 
The determination of q T at 0(M) gives access to the transverse momenta. Experimentally one component of q T is 
found as the non-collinearity of the produced particles K\ and K 2 in the plane perpendicular to the colliding particles 
Pi and P 2 , outlined in detail in Ref. Q. 

Accessing intrinsic transverse momenta in most cases requires a careful study of azimuthal dependence in high 
energy processes. Although the effects are in principle not suppressed by powers of the hard scale in comparison 
with the leading collinear treatment, it requires measuring hadronic scale quantities (transverse momenta) in a high 
momentum environment. Symmetries, in particular time reversal (T) invariance play an important role: 

• The theory of QCD is T- invariant. This makes it sensible to distinguish quantities and observables according 
to their T-behavior. 

• For distribution correlators involving plane wave hadronic states in the definition, combination of the T-operation 
and hermiticity, shows that the collinear correlators $(x) and T(x) must be T-even. For the TMD correlators, 
however, the T-operation interchanges & + \x,p T ) «-> &~}(x,p T ) (and similar relations for gluon TMD correla- 
tors). This allows to construct T-even and T-odd combinations. 

• For fragmentation functions the appearance of an hadronic out-state in the definition, prohibits the use of T- 
symmetry as a constraint and one has always both T-even and T-odd parts in the correlator (one can refer to 
T-even or T-odd in as far as the operator structure is concerned, referred to as naive T-even or naive T-odd). 

• In a scattering process, in which T-symmetry can be used as a constraint, single spin asymmetries would be 
forbidden. In fact the only real example of this is DIS (omitting electromagnetic interaction effects). For hadron- 
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FIG. 3: The paths in the gauge links in Q q (x,p T ). They involve a loop C/' D ' = U^~^ , which in the path shown in the left 
figure is closed (color-trace) and in the right figure is followed by a [+]-path. We will use the shorthand notations $[( D ) + ' and 
<E>' D+ ' respectively with (□) indicating the color-tracing and averaging. 



hadron scattering, e.g. the Drell-Yan process, one has a two-hadron initial state and only the assumption of a 
factorized description would imply absence of single spin asymmetries. We now know that this assumption is not 
valid, even not at leading order! Similarly for processes with identified hadrons in the final state T-invariance 
does not give constraints. 

• At leading order in a s , however, it is possible to connect single spin asymmetries (T-odd observable) to the 
T-odd soft parts, since the hard process will be T-even at this leading order. Collins and Sivers effects as 
explanation for single spin asymmetries are the best known examples. 



IV. THE TMD MASTER FORMULA 



The description of a hard process is obtained by writing down hard processes involving quarks and gluons and 
connecting these to the soft parts corresponding to initial state hadrons and observed hadrons in the final state. In 
the region where the hadrons are separated far enough in phase space (Pi ■ Pj ~ Q 2 , as discussed in the Introduction) 
one can have a soft part for each of the hadrons. For the determination of the relevant multi-parton matrix elements 
that need to be included in the calculation one can use the twist analysis alluded to in Section [TT] At leading order one 
needs the leading twist matrix elements (iptp) and (A T A T ), but also the multi-parton matrix elements (ip A + . . . A + ip) 
and (A T A + . . . A + A T ) (all having the same twist). The various matrix elements are resummed into color-gauge 
invariants combinations. Inclusion of the transverse pieces at infinity requires a careful analysis Q. The links that 
arise are process-dependent. They arise from diagrammatic contributions where collinear gluons A ■ n belonging to 
a particular soft part are attached to parton lines belonging to a different soft part (which are precisely the external 
parton lines of the linking hard subprocess). The link structure, thus, is not affected by inclusion of QCD corrections. 
On the other hand, the link structure depends on the color-flow in the specific diagram. 

The resulting expression for a hard cross section at measured q T is 

£ <e (B %i,^)*p ( ^ (i7) 

D.abc... 

where the sum D runs over diagrams distinguishing also the particular color flow and abc ... is the summation over 
quark and antiquark flavors and gluons. All Dirac and Lorentz indices, traces etc. are suppressed. The ellipsis at the 
end indicate contributions of other hard processes. 

We illustrate this master formula in the following example, taken from Ref. ||, describing the contribution in a 
hard scattering process coming from the qq — > qq subprocess (with both quarks having the same flavor). There are 
four diagrammatic contributions (see Table fl}, the result of which can be denoted as o q ^U qq with D running over the 
diagrams. For the first diagram, there are two different possibilities for the color flow, which absorbing the overall 
color factor in a have strengths (N 2 + 1)/(N 2 — 1) and —2/(N 2 — I), respectively. The second diagram has the same 
color flow possibilities. The third and fourth diagrams also have identical color flow possibilities but different from 
the first two diagrams. In this case each of the diagrams contributes two terms to the sum in Eq. II7[ e.g. the first 
diagram yields 

1,0 $[(□)+] (i) $[(□)+] ( 2 ) N c+ 1 & [D 1 ] A [C)-t] (1 / )A [C)-t ](2 / ) 
q y^/^g nj2 _ i qq-*qq q v > q \ > 



d 2 q T q w 9 w N 2 - 1 



+ <lf +1(1) *P+I(2) a%i qq A^-V) Af-V) + • • • , (18) 



where the underbraced items are separate a ^ -entries in the I?-summation of Eg. 1 1 71 
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oc ( 7/(0) { ™L - £f±I } m ) 



TABLE I: Gauge-links appearing in the soft parts connected to the qq — > subprocess depend on the specific diagrams. The 
paths for & q (x,p T ) are shown in Kig. [3] 



INTEGRATED AND WEIGHTED CROSS SECTION 



The results for cross sections after integration over the transverse momenta q T involve the path-independent inte- 
grated correlators Q(x) rather than the path-dependent TMD correlators & c ( D »(x,p T ). Thus, from Eq. [T7]one gets 
the well-known result 



E 

abc. . . 



where 



$aOi) ^b(x 2 ) a ab ^ c ..A c ( Zl ) 



/ j ao->c... 



(19) 



(20) 



is the partonic cross section. 

Constructing a weighted cross section (azimuthal asymmetry) by including a weight q" in the g T -integration leads, 
with the help of the relation between the observable q T and the intrinsic transverse momenta (e.g. the relation 
q T = p T — k T in SIDIS), to soft correlators of the form 

*Z lC] (x)= J d 2 p T P ^ [C] (x,p T ). (21) 

These still contain a path dependence, so Eq. [T7] cannot be simplified immediately. However, it turns out that the 
correlator in Eq. [5l] can be expressed as 



(22) 



Here $a(x) is a collinear correlator containing matrix elements with T-even operators, while ^>q(x,xi) is a collinear 
correlator with a structure like the quark-gluon-quark correlator shown in Eq. U involving the gluon field F na . In 
Eq. [22] one needs the zero-momentum (x\ — 0) limit for the gluon momentum. This matrix element is known as the 
gluonic pole matrix element. The operators involved are T-odd. Both collinear correlators on the RHS in Eq . |2"21 are 
link-independent. The gluonic pole factors Cq multiplying the gluonic pole correlator in Eq. I22[ however, do depend 
on the gauge link. They can be easily calculated. We have for instance = ±1, C^ + ' = 3 and C ( K D ' )+ ' = 1. Thus, 
one can write for the single-weighted cross section 

d 2 a 



d qr q™ 



d 2 qr 



D,abc... 



E 

abc. . 



$g a (xi) $ 6 (lE2)<7ah^c..A c (zi) ... + ... 

T^Ga^^l) ^b(x 2 ) 0-[o]6->c...Ac(^l) • ■ ■ + 



E 

abc... 



(23) 
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where the first term is multiplied by the normal parton cross section (Eq. I20|) and the second one involves the gluonic 
pole cross section, 



~ , _sr r iu(c(D))]^[D] (9A] 

0[a]b^c... ~ 2^°G a ab~+c... ■ \ l V 

D 



Noteworthy is the fact that these gluonic pole cross sections like the normal partonic cross sections also constitute 
gauge invariant combinations of the squared amplitudes. While for the electroweak processes like SIDIS and DY one 
has a simple factor, (T^ q ^i q = +ae q ^eq and cF[ q ]q^d = —Vqq^th the result for qq — > qq is more complex, 

= + (25) 

*w™ = + ^ 2] ) - frf (- [D3] + * [Di] ) . (26) 

where ct' 15 ^ refer to the contributions coming from the diagrams in Table |TJ Actually the results simplify in the limit 
N c — > oo in which case the color flow is unique for each diagram. Explicit results for gluonic pole cross sections are 
given in Ref. Q. 

The approach to understand T-odd observables like single spin asymmetries via the TMD correlators and the non- 
trivial gauge link structure unifies a number of approaches to understand such observables, in particular the collincar 
approach of Qiu and Sterman [|| and the inclusion of soft gluon interactions by Brodsky and collaborators [3] ■ Although 
the treatment of fragmentation correlators also separates into parts with T-even and T-odd operator structure, gluonic 
pole contributions (T-odd parts) in the case of fragmentation might vanish. Indications come from the soft-gluon 
approach Q and a recent spectral analysis in a spectator model approach [Iol |. 

VI. UNIVERSALITY 

Clearly the TMD master formula in Eq.[T7] breaks universality in the sense that one needs to know TMD correlators 
& u \x,p T ) with all sorts of gauge links U. Certainly it would be desirable to perform a study of the effects on the 
soft correlators caused by more complex gauge links than the simple ones given in Fig. O 

A useful procedure is to rewrite the TMD correlators in terms of T-even and T-odd correlators constructed from 
those in Fig. [5] and a residual or junk part, 

m {x,p T ) = i$[ cvcn l + \ C l " ] $[° dd l (x,p T ) + ( X) pr) ) (27 ) 

which by construction leads to 

&l«™4(x)=$(x), $^ [ovcn] (x) = $g(a;) (28) 

$ [odd] (x) = 0, $g [odd] (x) = vr x) (29) 

6<S> [u] (x) = Q, 8®% [V] {x)=Q. (30) 

For the simple quark correlators one has = 5$H = and the T -even and T-odd combinations are 

& e ™\x,p T ) = | (*M(x )Pt ) + 4> H (z,p T )) (31) 

&° d $(x,p T ) = \ ($M(x )Pt ) - t> [ - ] (x,p T )) . (32) 

Then one finds for instance 

& n +\x,p T ) = ±& even l(x,p T ) + §& odd ^x,p T )+5& D +^x,p T ), 

= 2<& [+] {x,p T ) - $ [ - ] {x,p T ) + 6$ [a+] (x,p T ). 

In Ref. Q it was noted that some further non-trivial simplification occurs for quark junk TMD while also the simple 
gluon correlators (Fig. ^ can be regrouped into T-even and T-odd combinations, 

T^ cn \x,p T ) = irl+'+Wx) + \V^-\x,p T ), (33) 
T [ ° dd] (x,p T ) = \ T^(x,p T ) - \ T^-\x,p T ), (34) 

rt° dd] (x, PT ) = ±r^-\ x ,p T ) - ^-^(x, Pt ). (35) 
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The two T-odd correlators reduce to the two three-gluon gluonic pole correlators Tq(x,x) and ^(1,1), which differ 
in the way the three color octets are coupled to a color singlet. 

The contributions S$ and ST make the non-universality explicit, which is a first step if one wants to study and 
possibly prove factorization in the case of TMD correlators. For phenomenological studies a reasonable first step is 
to omit the junk TMD, knowing that they will average to zero in a weighted asymmetry. This approximation can be 
applied immediately in the TMD master formula (Eq. [T7|) in which the color flow possibilities are dis ting uished. This 
master formula remains the basic starting point, used in recent analyses of photon-jet and jet-jet [l2j. production 
in hadron-hadron scattering. In the case of a linear weighting with the transverse momentum one can conveniently 
cast the result into folding of T-even and T-odd functions with normal and gluonic pole partonic cross sections, 
respectively. The procedure has been investigated for several processes (T3 j. in particular comparing the effect of 
using normal versus gluonic pole cross sections. 
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